Bayesian Algorithm for Unsupervised Unmixing of Hyperspectral Images Using a Post-nonlinear Model

This paper presents a nonlinear mixing model for hyperspectral image unmixing. The proposed model assumes that the pixel reflectances are post-nonlinear functions of unknown pure spectral components contaminated by an additive white Gaussian noise. The nonlinear effects are approximated by a polynomial leading to a polynomial post-nonlinear mixing model. A Bayesian algorithm is proposed to estimate the parameters involved in the model yielding an unsupervised nonlinear unmixing algorithm. Due to the large number of parameters to be estimated, an efficient constrained Hamiltonian Monte Carlo algorithm is investigated. The performance of the unmixing strategy is finally evaluated on synthetic data.